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A relativistic complex scalar boson field at finite temperature T is examined below its 
critical Bose-Einstein condensation temperature. It is shown that at the same T the state 
with antibosons has higher entropy, lower Helmholtz free energy and higher pressure than 
the state without antibosons — but the same Gibbs free energy as it should. This implies that 
the configuration without antibosons is metastable. Results are generalized for arbitrary d 
spatial dimensions. 
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INTRODUCTION 



In early works [l-3] on the relativistic ideal boson gas (RIBG) explicit Bose-Einstein conden- 
sation (BEC) critical transition temperature T r -formulae were derived for both the nonrelativistic 
and ultrarelativistic limits and specific-heat anomalies at T c were studied. In addition, Refs.|2, 3] 
considered all space dimensions d > and delved into the relation between d and various critical 
exponents. At sufficiently high temperatures, however, boson-antiboson pair production becomes 
appreciable and this was not accounted for. The first reports to include both bosons and anti- 
bosons appear to be Refs.jJ, Q| where high-temperature expansions for the various thermodynamic 
functions (pressure, particle-number density, entropy, specific heats, etc.) were derived. Extensive 
numerical work in d dimensions that does not rely on such high-temperature expansions was re- 
ported in Refs.y, Q|. In the elegant treatment of Ref.[8] with inverse Mellin transforms the specific 
heat anomaly of the RIBG at its BEC T c was found to be washed out when pair-production was 
included. The relationship between the BEC of the RIBG and spontaneous-symmetry breaking 
was explored in Refs.jj], j|; see also the rather complete Ref.ficj]. esp. §2.4. 

BECs are also of interest in cosmological and astrophysical contexts. In fact, increasing attention 
has recently been paid cosmological models that describe dark matter (DM) as a condensate phase 



of some scalar boson field 
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24j. Such models are competitive with the A cold dark matter 



(ACDM) model 25f] to explain observational properties of DM at cosmological and astrophysical 
levels. In particular, a scalar boson field with an extremely small mass of about 10 -22 eV can 



explain the cosmological evolution of the universe 
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-241], the rotation curves of galaxies [26], the 
central-density profile of low-surface-brightness galaxies 27 K the size of galactic halos [281 ]. and 
the amount of substructures in the universe [29(]. In Ref. [301] it was shown that a complex and 
self-interacting scalar boson field with a more realistic mass of about 1 eV in a BEC is also a viable 
DM candidate. Moreover, in this model no fine tuning of the scalar-field energy density at early 
times is required and the condensate formation is due to self-interactions. 

Indeed, BECs are of interest in the context of quantum gravity. In Refs. 



42 



431 ] some of 



us have shown that Planck-scale deformations of the energy-momentum relation that natural 



emerges in many quantum-gravity theories (for an excellent nontechnical overview see Ref. 441 ]) 



may affect the properties of low-temperature BECs. In particular it was shown that a Planck-scale 



induced deformation of the Minkowski energy- momentum dispersion relation E ~ ^/m 2 c 4 + p 2 c 2 + 
£ t mcp/2M p , where m is the mass of the bosons, M p the Planck mass and £ a dimensionless 
parameter, produces a shift in the condensation temperature T c of about AT C /T C ° ~ 10 -6 £i in 



3j, 7 3 U [3|, ?3 Na [33J, i H [3|, |5Rb [3j, |He [36j, gK |37|, i 3 Cs 



typical BECs such as %lRb 

n n 

3a |. and llCr [39(]. The quantum gravity induced shift in T c makes possible to upper-bound the 



deformation parameter as |£| < 10 4 with recent ultra-precise measurements of T c as, e.g., in fgK 
j4o| . In Ref.[4^. 13] it is also discussed how to enlarge AT C /T C ° thus improving the bound on £ 
and hence realize an ad hoc experiment accomplish this. Finally, the Planck-scale induced shift 
in T c is compared with similar effects due to interboson interactions and finite-size effects. These 
results open a new possibility for a quantum gravity phenomenology based on low-temperature 
condensates, so that BECs truly appear to be a frontier interdisciplinary research field open to 
many applications. 



We also stress how the effect of interactions as well as of finite-size effects might have observable 
effects on laboratory BE condensates. For example, in Ref.j^J the effect of interactions has been 
observed in fgK and a shift in T c measured as a function of the interboson s-wave scattering length 
a and data have been fitted with the second-order polynomial AT C /T® ~ b\ (a/Ay) + ^(ci/At) 2 
with b\ = —3.5 ± 0.3 and 62 = 46 ± 5, the second term being due to beyond-mean-field effects. 
However, in what follows we do not consider interaction nor finite-size effects as we focus on the 
ideal Bose gas. Such effects are being investigated. 



Here we study the metastability of a BEC that does not contain antibosons. A motivation 
is given in §11. In Ref. 4l| the properties a RIBG in terms of the Helmholtz free energy with 



antibosons included was discussed and shown to be a state with a lower Helmholtz free energy 
than that without antibosons. In §111 we generalize this result by comparing two different BECs, 
with and without antibosons, but with the same total number of particles and at the same finite 
temperature. Such states are related by a thermodynamic transformation ensuring that they are 
meaningfully comparable. In particular, we rely on the law of nondecreasing entropy for isolated 
systems. In §IV we conclude that the state with antibosons has greater entropy and lower Helmholtz 
free energy and which is therefore the stable state, while the state without antibosons is metastable. 
In §V we derive the expression of the pressure of the BEC in equilibrium with a thermalized gas 
of bosons and show that the state with antibosons has higher pressure. As an overall check we 
calculate the Gibbs potential in both cases and show that it is the same, as expected. Lastly, in 
§VI we generalize results for arbitrary d > dimensions, integer or not. We conclude in §VII. 
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II. MOTIVATION 

To study the relative stability of the two states with and without antibosons one should compare 
their entropies. For a meaningful comparison the two states must be at the same temperature, 
volume and number density. This is guaranteed in what follows. Consider a system composed of 
two heat reservoirs R\ and R2 at temperatures T\ and T2, respectively, with T\ <C I2, and a gas of 
iV bosons B of mass rest mass m contained in a volume V with a number-density n = N/V . The 
reservoirs are much larger in volume than the boson volume V so they can be placed in thermal 
contact with the boson gas without appreciably changing its temperature, or they can be isolated 
from the boson gas. Assume that T x < Tf < T^ S , where Tf and Tf B are the boson gas BEC 
critical temperatures without and with antibosons B, respectively, so that at T\ the boson gas 
itself is a BEC. Assume also that ksT2 » mc 2 , with c the velocity of light, and that ksT\ <C mc 2 . 
Initially, the gas contains only bosons and is in thermal equilibrium with the first reservoir R\ at 
temperature T x . Since k^Ti <C mc 2 at this temperature any antibosons present are negligible. The 
boson gas is then isolated from the reservoir R\ and placed in thermal contact with the reservoir 
i?2- After awhile the boson gas reaches thermal equilibrium at temperature Ti- Since fce^ ^> mc 2 
antibosons are created substantially by pair-production so that the equilibrium state now also 
contains antibosons B. Finally, the boson gas is isolated from the reservoir R2 and placed in 
thermal equilibrium with the reservoir R\ so that the final temperature of the boson gas is T\. 
The question concerning the metastability of the state without antibosons can be formulated in 
the following way: at the end of the process just described does the boson gas contain antibosons 
or does it go back to the initial state without antibosons? To answer this one must calculate the 
entropy variation A5/ ot of the whole system (boson gas and reservoirs) a final state of the boson 
gas with antibosons, and compare it with the entropy AS^ of a final state without antibosons. 
This question is addressed and resolved in §IV where we show that ASl ot > AS^ t so that the state 
without antibosons is metastable. 

We first calculate the main thermodynamic functions in both cases, with and without anti- 
bosons. 

III. ENERGY DENSITY AND HELMHOLTZ FREE ENERGY BELOW BEC T c 

We consider two gas systems, one with only bosons B and a second one containing also anti- 
bosons B. They are both as the same temperature T < < T C BB , where Tf and T C BB are the 
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condensation temperatures of these two systems without and with antibosons, respectively. We 
first write down explicit expressions for internal energies and number densities and then proceed 
to calculate their Helmholtz free energy. 

Since T < T B < T BB the condensate forms in both a system containing only bosons B as 
well as in a system containing also antibosons B. At such temperatures T the chemical potential 
fi ~ mc 2 in a RIBG whose energy i?(p)-momentum p dispersion is E(p) = \Jp 2 c 2 + m 2 c A . For a 
gas containing only bosons the number density is 

n = n + (h^ir 2 )" 1 / p 2 dp -— — ^ 7 (!) 

Jo+ ex P \P( E (P) - mc 2 )} - 1 

where f3 = 1/ksT. The net internal energy per unit volume V is 



U B (n,T,V) 2 ,,.„„ 2N , f°° 2 , E(p) 

J + exp \fi(E(p) - mc 2 )] - 1 



V 
Here 



= 1 1 (O) 

U ° V exp[/3(mc 2 - //)] - l' 1 j 



Combining these equations leaves 



U B (n,T,V) 2 „o 2x i Z" 00 2 , %)-mc' 

- mc 2 n + n 3 27r 2 - 1 / P 2 dp f — ^ f^— r- 4 

Jo+ ex P \p{E(p) - mc 2 )} - 1 



V 

When antibosons are included the number density n is 



n 



n + (hV)- 1 1 



p dp 



0+ Lexp \P(E(p) - mc 2 )] - 1 exp [/3(E(p) + mc 2 )] - 1 



(5) 



so that 



U (n, T, F) 



/"OO 

= mc 2 n +(?i 3 27r 2 )" 1 / p 2 dpE(p) 



+ 



exp — mc 2 )] — 1 exp \j3(E(p) + mc 2 )] — 1 

(6) 



Combining these two equations gives 



U 



BB 



( n ' T ' V l=mt?n+(hV)- 1 ^ p 2 dp 



V 



E( P ) 



m c 



exp [p(E(p) 



mc^ 



+ 



E( P ) 



m c 



exp [(3(E(p) + mc 2 )] - 1 

(7) 
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The Helmholtz free energy per unit volume without antibosons is then 

F B (T,V,n)/V = mc 2 n + k B T(h 3 27r 2 )- 1 p 2 dp In [l - exp [0 [mc 2 - E(p)] )] . (8) 
In the case with antibosons one has 

poo 

F BB (T,V,n)/V = mc 2 n + k B T{h 3 2ir 2 )- 1 / p 2 dp{ln [l - exp [0 (mc 2 - E(p))]] + 



ln[l-exp[-/3(mc 2 + £(p))]]}. (9) 
From ([HI)-® it also follows that 

F BB (T,V,n) - F B (T,V,n) = Vj^^ J p 2 dpln [l - exp (-/3 [mc 2 + E(p)] )] < 0. (10) 
Therefore the state containing antibosons has a lower Helmholtz free energy. This same result 



4l| except that here the Helmholtz free energies are compared at the same 



was found in Ref. 
temperature T. 



IV. ENTROPY 

Here we calculate the entropy of the boson field with and without antibosons. This is then 
used to determine the entropy variation in the thermodynamic transformation described in ^Hlto 
conclude that the state containing antibosons is the stable state while the state without antibosons 
is only metastable. If only B bosons are considered, the entropy follows from 

TS B (T,V,n) = U B (T,V,n) - F B (T,V,n) (11) 
where the internal energy per unit volume is given by ([!]). Whence 

S°P, V,n)IV . k B (HV)- 1 f^M g^Hl^p^ -ln [1 - exp (P - JSfp)])] }. 

(12) 

If antibosons are included the entropy follows from 

TS BB (T, V, n) = U BB (T, V, n) - F BB (T, V, n) (13) 
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where the Helmholtz free energy is given by ((9|). Using the latter and ([7]) one gets 

q BB fTV wv , (h s 9 2,-1 r 2, r (3(E(p)-mc 2 ) f3(E(p) + mc 2 ) 

S (T, V, n) V = k B {n 2tt ) / p dp{ ^ — - + =r= — -+ 

j + exp [p(i?(p) — mcr)J — 1 exp [p{E(p) + mc z )\ — 1 



In [l — exp [P {mc 2 - E(p))]] - In [1 - exp [-/3 (mc 2 + S(p))]]}. (14) 



One can nowcompare the entropies of the two states with and without antibosons. From (|12p and 
([HI) one easily finds that 



poo 

S BB (T,V,n)-S B (T,V,n) = ^(fr^yr 2 )- 1 / p 2 dp{fi (E(p) + mc 2 ) /[exp [/3(£(p) + mc 2 )] -1] + 



— In [l — exp (-0 [mc 2 + E{p)} )} } > (15) 

so that the state without antibosons being less entropic is thus metastable. 

Now consider the thermodynamic transformation described in £jll]and calculate the total entropy 
variation of the boson field plus that of the two reservoirs. This enables one to decide if the final 
state will contain or not antibosons. The final state of the whole system (boson field plus reservoirs) 
turns out to be more entropic one which in turn implies that the state with antibosons is the stable 
state while the state without antibosons is only metastable. If the final state of the gas also contains 
antibosons, the entropy variation of the gas in the thermodynamic transformation described in 21] 
is 



AS I gas = S BB (T 1 ,V,n)-8 B (T 1 ,V,n)>0 (16) 
(which is the same as ()15p ). while the entropy variation of the two reservoirs is 



AS* = = U BB {T2,V,n)-U B {T u V,n) 

l T Ti { ) 



, AQ 1 = U^T 1 ,V,n ) -U^ i T 2 ,V,„) 

2 T 2 T 2 V ; 

Hence, the total entropy variation is 
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ASL = AS^+ASi+ASi = A^ as +t/ BB (T 2 , V, n) ( - - - ) + 



gas 



1 1 \ , U BB (T u V,n) U»(li,V,n) 



B i 



T 2 



> 



> A5^ + 



c/ B5 (r 2 ,y,n)-^(r l5 y,n 



1 1 , 



(19) 



The transformation is thus allowed but is irreversible. 

The net entropy variation in the thermodynamic transformation of $11] when the final state is 
without antibosons is thus 



AS£L = 



(20) 



while the entropy variation of the two reservoirs is 



AS U = U BB (T 2 ,V,n)-U B (T u V,n) > Q 



(21) 



as ii = U B (T u V,n)-U BB (T 2 ,V,n) 
2 T 2 



(22) 



Hence, the total entropy variation is 



U BB (T 2 ,V,n) -U B (Tx,V,n) 



(23) 



Again, the transformation is allowed but is irreversible and its only effect is a heat transfer between 
the two reservoirs Ri and R 2 . 

We can now compare the two entropy variations AS[ ot and AS// 4 . One has 



AS[ ot - AS& 



ASg as + 



U BB {T l ,V,n)-U H {T 1 ,V,n 
T 2 



> 



(24) 



so that ASt ot > A Slo t , i.e., the entropy variation is greater in the case in which the final state 
contains antibosons. Therefore, the final state of the whole system of the gas plus the two reservoirs 
is more entropic when the gas contains antibosons in the final state. Again this means the state 
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without antibosons is metastable and that the final equilibrium state described in £JTT] is the one 
containing antibosons. 

We remark that if the boson field were a real scalar field it would not admit antibosons and 
in this instance the state containing only bosons B is the only possible one and therefore it is not 
metastable but rather a stable state. One thus concludes that if antibosons are allowed, namely 
if the scalar field is complex as assumed here, the state with antibosons is allowed and this state 
will be the stable one while the state without antibosons will only be metastable. 

V. PRESSURE 

Following the same procedure we introduce the pressure P as a function of n, V and T . 
Specifically, if no antibosons are present 

/•oo 

(3VP B = - In [1 - exp [/3(fx - mc 2 )]] - (h 3 2ir 2 )- 1 V / p 2 dpln [1 - exp - E(p))}} . (25) 
We rewrite ([I]) as 



1 exp[/3(^-mc 2 )] 

n = n-n + = — — ( 26 

VI — exp[p(^ — mc 1 )] 



where 



f°° 1 

n+ = (h 3 27r 2 )- 1 / p 2 dp r-TTTTT-, ^ (27) 

+ ' Jo+ exp \p(E(p) - mc 2 )] - 1 V ' 

is the number density of noncondensate (or excited) bosons. Below the condensation temperature 

^ ~ mc 2 so that exp [(3(fi — mc 2 )] ~ 1 apart from small corrections 0(1/V) which vanish in the 

thermodynamic limit V — > oo. We can write the logarithm in the rhs of ()25p as 



In [l — exp [f3(n 



mc 2 )]] 



In [V(n — n + )] 



(28) 



whence 



(3 VP 



In 



V n 



(n 3 2vr 2 ) 



p 2 dp 



exp \J3 (E(p) — mc 2 )] — 1 J 



V 



(h?2n 2 ) 



/•oo 

/ p dp In [1 — exp [f3 (mc 2 - E(p))]] . (29) 
./()• 
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Dividing through by V gives the first term on the rhs proportional to V 1 In V which also vanishes 
in the thermodynamic limit, so one gets 

p b = -T^iy £p 2d P in t 1 - ex p^ ( mc2 - ■ ( 3 °) 

When antibosons are included the pressure is given by the relation 

1 f°° o . ( 1 1 



f3VP BB = - In 



Vn — I p 2 dp 



h?2n 2 J + " r V ex P \P(E(p) ~ mc2 )\ ~ 1 ex P [P(E(p) + mc2 )\ ~ 1 



+ 



V 



{h 3 2n 2 ) 



/ p 2 dp [In (l - exp[/3 {mc 2 - E(p))}) + In (l - exp[-/3 (mc 2 + E(p))]) ] . (31) 



Since the first term in ()31[) is negligible in the thermodynamic limit, the final result is 

P BB = -T^y /Jp 2 ^ [1 - exp [p [mc 2 - E(p))}] + In [l - exp [-0 {mc 2 + E(p))]]) . 

(32) 

Comparing (f30l) with ([32]) it becomes evident that when antibosons are included the pressure is 
greater than in the case without antibosons. Indeed, one has 

P bb _ pB = - Jg^ £Vdpln [1 - exp [-p {mc 2 + E(p))]] > 0. (33) 



As a final overall check, comparison of equations (|10p and (|33p shows that the state with 
and without antibosons have the same Gibbs free energy G(P, T) = F + PV = fiN, namely 
G ■ (P, T) = G B {P,T), as must be the case since the net number of particles N and chemical 
potential fj, are the same. 

VI. GENERALIZATION TO d SPATIAL DIMENSIONS 

Here we generalize the thermodynamic potentials to arbitrary d > spatial dimensions, integer 
or not. The result is confirmed that the state with only bosons is metastable while the state 
with both bosons and antibosons is stable. To motivate this section we recall that spaces with 
dimensionality different from d = 3 are considered in many physical contexts, e.g., in quantum 
gravity (see Ref. 



441 ] for a review). In other areas, e.g., Mandelbrot (Ref. 



451 ]. p. 85) cites an 



empirical fractal dimension d = 1.23 for the distribution of galaxies in the observable universe. 
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We first calculate the thermodynamic functions with and without antibosons. Assuming a real 
nonnegative number d of spatial dimension, the sum over momentum now becomes 

<"> 

where Q d is the solid angle in d dimensions and the system volume is L d . 
If no antibosons are present the number density n is 

POO 

n = n + n d (2nH)- d / p d ~ x dp [exp \p (E(p) - fj,)\ - l]' 1 (35) 
Jo+ 

where 

n = [V (exp [0 (mc 2 - M )] - l)] _1 . (36) 

As before, no = No/L d is the number density of zero-momentum p = bosons within a d- 
dimensional volume V = L d . The internal energy per unit volume is 

U B (T,V,n)/V = nmc 2 + n d (^r d H p d ~ X dp f/gj (37) 

7 + exp [/3 (£(p) - m)] - 1 

The Helmholtz free energy is 

;»oo 

F B (T,V,n)/V = mc 2 n + k B Tn d (2irh)~ d p d ~ x dp In [1 - exp [/? (/i - E(p))]] (38) 

while the entropy per unit volume is now 

S B (T, V, n)/V = k B n d (2nh)- d H p d ~^dp{ P f ^ ~ ^ - In [1 - exp [/3 (u - E(p))]]} 

Jo+ exp [/3 (E(p) - (j.)] - 1 

(39) 

with the term (k B /V) In [l — exp [/? (/i — mc 2 )]] being negligible in the thermodynamic limit. 
If antibosons are present the number density is 



poo 

n = n + n d (2irh)- d / p^ 1 dp 



1 

+ 



_exp \p(E(p) - m)] - 1 exp [9(JS?(p) +/*)]- 1 
where no is still (|36|) . The internal energy per unit volume is 

U BB (n,T,V) _ 2 , n ,„ fc Hf J-uJ £(;p)-m C 2 , %)+mc 



(40) 



poo 

mc 2 n+n d {2nhy d / p d - l dp 



V 

The Helmholtz free energy per unit volume is 



+ 



exp \p(E(p) — mc 2 )] — 1 exp [(3(E(p) + mc 2 )} — 1 

(41) 



poo 

F BB (T,V,n)/V = mc 2 n + k B TQ d (2iThy d / p d_1 dp{ln [l - exp [j3 [mc 2 - E(p))]] + 
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In [1 -exp [-/3(mc 2 + E(p))]]}. (42) 



Finally, the entropy per unit volume becomes 

("OO 

/o P dp ^laq> \P(E(p) - mc 2 )} - 1 ' exp \J3(E(p) + mc 2 )} - 1 



B5 d /- 00 d _! /3(g(p)-mc 2 ) /3(£(p)+mc 2 ) 
6 (r,V,n)/K = fcBO d (27rn) / p dp{ ^ -H gvi -+ 



— In [l — exp [P {mc 2 - E(p))]] - In [1 - exp [-/3 (mc 2 + E(p))]]}. (43) 

At this point it is easy to generalize to arbitrary d the result that the state with only bosons 
is stable, just generalizing (f24"|) . For example, by use of (j37ti38tf39ti41ti42tf^3|) one easily generalizes 
(fT0|) and obtains 

F B *(T,V,n)-F B (T,V,n) = Vd ^^ £ p^ 1 dp In [l - e W [mc 2 + E(p)] )] < (44) 

for the difference of the Helmholtz potential. Again, for arbitrary d the state containing antibosons 
has a lower Helmholtz free energy. Then one can generalize (|15p and obtain 



S BB (T, V, n) - S B (T, V, n) = ^0 £ p d ~ l dp{P {E{p) + mc 2 ) /[exp [(3(E(p) + mc 2 )] - 1]+ 



- In [1 - exp [-p [mc 2 + E(p)])] } > (45) 

so that the state with antibosons is more entropic also for arbitrary d. 

Proceeding in the same way one verifies that the relations resumed in Eq.s ([15M23|) are still valid 
and therefore Eq. (j24p is also valid for arbitrary d. Therefore one concludes that, also for arbitrary 
d, the state without antibosons is metastable. 

VII. CONCLUSIONS 

The metastability of a Bose-Einstein condensate (BEC) that does not contain antibosons was 
studied for the relativistic ideal Bose gas (RIBG). In particular, the Helmholtz free energy with 
both bosons and antibosons was shown to be a state with a lower Helmholtz potential than that 
without antibosons. This was done with the same number of particles and at the same finite 
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temperature below the BEC critical temperature. Both states were found to be related by a ther- 
modynamic transformation ensuring that they are meaningfully comparable. In addition, relying 
on the principle of nondecreasing entropy for isolated systems we found that the state with an- 
tibosons has greater entropy and is therefore the stable state, while the state without antibosons 
is metastable. The pressure of both systems was calculated and found to be higher for the state 
with antibosons than for the state without them. We also confirm that the two states with and 
without antibosons have the same Gibbs free energy, as expected. Lastly, results were generalized 
for arbitrary dimensions d > 0, integer or not. 
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